Superfluidity and Interference Pattern of Ultracold Bosons in Optical Lattices 
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We present a study of the superfluid properties of atomic Bose gases in optical lattice potentials 
using the Bose-Hubbard model. To do this, we use a microscopic definition of the superfluid fraction 
based on the response of the system to a phase variation imposed by means of twisted boundary 
conditions. We compare the superfluid fraction to other physical quantities, i.e., the interference 
pattern after ballistic expansion, the quasi-momentum distribution, and number fluctuations. We 
have performed exact numerical calculations of all these quantities for small one-dimensional systems. 
We show that the superfluid fraction alone exhibits a clear signature of the Mott-insulator transition. 
Observables like the fringe visibility, which probe only ground state properties, do not provide direct 
information on superfluidity and the Mott-insulator transition. 
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Ultracold atomic gases in optical lattices provide a 
unique framework for the experimental study of funda- 
mental quantum phenomena in interacting many-body 
systems. This is especially true for the exploration 
of quantum phase transitions such as the superfluid to 
Mott-insulator transition observed in a recent pioneering 
experiment Q. The remarkable degree of experimental 
control over all the relevant parameters — density, inter- 
action strength, lattice geometry and dimensionality — 
allows much more detailed studies of the complicated 
mechanisms behind quantum phase transition than con- 
ventional solid state systems. 

It is clearly the case that the most important quan- 
tity for characterizing the superfluid-to-insulator transi- 
tion is the superfluid density or superfluid fraction /g. 
The aim of this paper is to set up the general theoretical 
framework for the calculation of the superfluid fraction 
within the Bose-Hubbard model and to compare it, on 
the formal level, with quantities being measured at the 
moment. These include the interference pattern after bal- 
listic expansion, the quasi-momentum distribution as well 
as the number fluctuations which are important in appli- 
cations. We perform exact numerical calculations for the 
Mott-insulator transition in an one-dimensional system 
to demonstrate the relationship, and more importantly 
the differences, between the superfluid fraction and vari- 
ous ground state observables, most notably the visibility 
of the interference pattern. 

Superfluidity. The concept of superfluidity is closely 
related to the existence of a condensate in the interacting 
many-boson system Formally, the one-body density 
matrix p^^\x,x') has to have exactly one macroscopic 
eigenvalue which defines the number of particles A^o in 
the condensate; the corresponding eigenvector describes 
the condensate wave function 0o(^) = e'^(=^)|0o A 
spatially varying condensate phase 9{x) is associated with 
a velocity field for the condensate by 

z/o(f) = - V0(f) . (1) 
m 

This irrotational velocity field is identified with velocity 
of the superfluid flow Vs{x) = vq{x) j2j and enables us 



to derive an expression for the superfluid fraction /g = 
Ng/N . Consider a system with a finite linear dimension L 
in ei-dircction and a ground state energy Eq calculated 
with periodic boundary conditions. Now we impose a 
linear phase variation 9(x) — Qxi/L with a total twist 
angle 8 over the length of the system in the ei-direction. 
Technically, this can be achieved by introducing twisted 
boundary conditions of the form ^'(xi, Xi + Lei, ■■■) = 
giO \i)(^xi, ...^Xi, ...) with respect to all coordinates of the 
many-body wave function. The resulting ground state 
energy Eq will depend on the phase twist. For very small 
twist angles Q <^ tt the energy difference Eq — i?o can 
be attributed to the kinetic energy Tg of the superflow 
generated by the phase gradient. Thus 

- ^0 - - \mNf,,vl , (2) 

where mNfs is the total mass of the superfluid compo- 
nent. Replacing the superfluid velocity Vs with the phase 
gradient according to Eq. leads to a fundamental 
relation for the superfluid fraction Q 
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Hence the superfluid fraction can be interpreted as a mea- 
sure for the stiffness of the system under an imposed 
phase variation. This demonstrates that superfluidity 
is not a static ground state property but rather the re- 
sponse of the system to an external perturbation. We 
should note that this definition of superfluidity does not 
tell anything about the stability of the superfluid flow at 
finite velocities. 

We now transfer these findings to a lattice system de- 
scribed in the framework of the Bose-Hubbard model 
^, ^. For simplicity we restrict ourselves to a regular 
one-dimensional lattice composed of / sites described by 
the Bose-Hubbard Hamiltonian 

^ V ^ 

Ho = - J^(aTa,+i + h.a.) + i^Yl " 1) ' 

i=l i=l 

where a| creates a boson in the lowest Wannier state lo- 
calized at site i and rij — aj^a^. The first term describes 
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the hopping or tunnehng between adjacent sites with a 
tunnehng strength J, the second term characterizes the 
on-site two-body interaction with a strength ||, 0] ■ The 
hopping between the first and the last site of the lattice is 
included (7+1=1), which corresponds to periodic bound- 
ary conditions. The generalization to three-dimensional 
lattices is straight forward. 

In order to compute the energy of the system with an 
imposed phase twist G we map the twisted boundary con- 
ditions by means of a local unitary transformation onto 
the Hamiltonian ||, ^. This leads to a twisted Hamilto- 
nian of the form 
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with a modified hopping term that contains the so-called 
Peierls phase factors e^'®/^ The energy Eq in the 
presence of the phase twist is given by the lowest eigen- 
value of the twisted Hamiltonian using periodic bound- 
ary conditions. From the difference of the ground state 
energies Eq — Eq we obtain the superfluid fraction [nOl 
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now expressed in terms of the parameters of the Bosc- 
Hubbard model. 

One can get a more detailed insight into the depen- 
dency of the superfluid fraction on the structure of the 
eigenstates of the system by considering a perturbative 
calculation of the energy difference Eq — Eg . We expand 
the the twisted Hamiltonian up to second order in the 
twist angle Q thus, 
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Here we defined a current operator J = it/X]i(^i-i-i^j 
h.a.) and the usual kinetic energy or hopping operator 

T = —JJ2ii^l+i^i + h.a.). We can calculate the energy 
shift Eq — Eq caused by the perturbation Hport in second 
order perturbation theory. Retaining the terms up to the 
quadratic order in the twist angle O we obtain for the 
superfluid fraction using Eq. (o) 
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where the {v = 0,1...) are the eigenstates of the 

non- twisted Bose-Hubbard Hamiltonian Hq . The ground 
state expectation value of the hopping operator describes 
the first order contribution /s'^^. The sum over the ex- 
cited states involving the current operator constitutes the 
second order term fj^'^\ It is this second order term 
which introduces a significant dependence of the super- 
fiuid fraction on the excitation spectrum and thus goes 
beyond the static ground state properties of the system. 




FIG. 1: Superfluid fraction /s as function of the interaction 
strength V/J for filling N/I — 1 and different lattice sizes: 
/ = 10 (full line), 8 (dashed), 6 (dotted). The thin lines show 
the first order term /i^'. The vertical line marks the critical 
interaction strength (V/J)crit ~ 4.65 for the infinite system. 



Equation (|) corresponds to the Drude weight used to 
characterize the DC conductivity of charged fermionic 
systems jll]. This demonstrates that the phase factors 
appearing in the twisted Hamiltonian (^) can actually be 
realized experimentally, i.e., by an external electric field 
for charged particles or by some linear external potential 
or even by accelerating the lattice. 

To illustrate the generic behavior of the superfluid frac- 
tion as function of the interaction parameter V/ J, and 
hence the appearance of the Mott insulator phase, we 
solve the eigenvalue problems of the non-twisted Hamil- 
tonian (^) and the twisted Hamiltonian ^ numerically. 
We construct the corresponding Hamilton matrices us- 
ing a complete basis of Fock states |rti, ...jni) with all 
compositions of the occupation numbers . The ground 
state energies Eq and -Ee=o.i are obtained with an effi- 
cient iterative Lanczos algorithm The perturbative 
expression (H) is then used to separate the contributions 
of ground and excited states. 

Figure |l| shows the superfluid fraction fs for one- 
dimensional lattices with up to /= 10 sites and fixed 
filling factor N/I = 1. The superfluid fraction is 1 for 
the noninteracting system and decreases slowly for small 
V/J. In the region of the Mott transition fs drops rapidly 
and goes to zero in the Mott-insulator phase. The se- 
quence of curves for increasing system size shows a mod- 
erate size dependence around the onset of the insulator 
phase. One can extrapolate the curves to / ^ oo and 
finds a vanishing superfiuid fraction above a critical in- 
teraction strength which is in good agreement with the 
value (V/ J)ci-it ~ 4.65 obtained by strong coupling ex- 
pansion [g2| and Monte Carlo methods [Q. The thin 
lines in Fig. |^ depict the isolated first order contribution 
fj^^ = -2]^(^'o|T|«'o) to the superfluid fraction (|), 
which is just the reduced expectation value of the hop- 
ping operator. This quantity decreases much slower than 
the total superfluid fraction. Thus even deep in the Mott 
regime the hopping operator has a considerable expecta- 
tion value (up to 30% of its value in the non-interacting 
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system) although the system is already an insulator, i.e., 
the superfluid fraction is zero. The rapid decrease of the 
total superfluid fraction is largely due to the second or- 
der contribution /s*-^-*, which vanishes for small V/ J and 
shows a threshold-like behavior around the Mott transi- 
tion. The vanishing of the superfluid fraction in the insu- 
lating phase is generated by a strong cancellation between 
the first and second order term. This emphasizes that the 
superfluid fraction depends crucially on the properties of 
the excited states. 

Interference Pattern. The standard experimental ap- 
proach to investigate the state of the Bose gas in an opti- 
cal lattice relies on the matter- wave interference pattern 
after the gas was released from the lattice. How much 
can the presence or absence of interference fringes tell 
about superfluidity? 

In the simplest model of the expansion for the system 
after release from the lattice the effects of interactions 
are neglected. We can write the intensity observed after 
some time-of-fiight r at a point y as 



X(y) = (M'o|At(y)A(y)|M'o> 



(9) 



We assume that the Wannier functions w{x — ^i) can be 
described by Gaussians. The amplitude operator is given 
by A(j7) = -^Yli^iXiiy) a^, where xdv) denotes the 
Gaussian wave packet associated with site i after a free 
evolution for a time r. Since we are interested only in the 
generic features of the interference pattern we discard all 
terms related to the spatial structure of the envelope of 
Xi{y) a-nd only retain the phase terms. This leads to 



(10) 



where (l)i{y) is the total phase acquired on the path from 
site i to the observation point y. In the far-field ap- 
proximation we can assume a constant phase difference 
S(j){y) = 4>i+i{y) — <t'i{y) for adjacent sites (gravity ne- 
glected). Calculating the intensity (^ as function of the 
phase difference 50, using the far-field form of the am- 
plitude operator (p^, leads to the following expression: 
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N + '^Bd cos{d6<j)) 



In the last step we rearranged the double summation into 
a sum over the hopping distance d — j — i. The coefS- 
cients Bd are given by the expectation values of the dth 
neighbour hopping operators 



Bd = Y,{'^o\4+d^^+4^^^ 
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FIG. 2: Intensity T{Scj>) as function of the phase difference 
54> for a system with 7 = = 10 and V/J = (full line), 5 
(dashed line), 10 (dotted line). 




FIG. 3: Visibility V of the interference fringes (full line) and 
number fluctuations a (dashed) as function of the interaction 
strength V/J for a system with J = A'^ = 10. 



The reader should note that the leading coefficient Bi is 
related to the first order term of to the superfluid frac- 
tion d) through Bi = 2{I - d) fj'^\ Clearly, there is 
no contribution corresponding to the important second 
order term /s*-^^ of the superfluid fraction, because the 
intensity (|^) involves only the ground state. Thus the in- 
terference pattern cannot provide full information on the 
superfluid properties, as it only measures the flrst order 
term of the superfluid fraction. 

Figure || shows the intensities I{S(t)) resulting from 
the exact numerical calculation for different interaction 
strengths. The interference peaks around (50 = 0, ±27r, ... 
correspond to the prominent peaks observed experimen- 
tally. Recall that terms describing the overall envelope 
were neglected in (^0|). With increasing V/ J the intensity 
2max of these principal peaks reduces. At the same time 
an incoherent background emerges such that the mini- 
mum intensity Imin between the principal peaks grows. 
Thus the interference fringes are increasingly washed-out 
and eventually only a flat intensity distribution remains. 

As a quantitative measure for the vanishing of the in- 
terference pattern the full line in Fig. ^ shows the fringe 
visibility V = (Xmax - 2:min)/(2'max + ^min) as fuuctlon 
of V/ J. In addition, the dashed curve shows the on-site 



4 



number fluctuations (Ti = ((n|) — (ni)^)^/^ of the ground 
state. Obviously, the visibility of the fringes has no im- 
mediate relation to the number fluctuations. For small 
interaction strengths V^/ J ^ 5 the visibility remains al- 
most constant at V « 1 whereas the number fluctuations 
drop to 0.5 in the same interval. 

A second observation concerns the relation with the 
superfluid fraction shown in Fig. [l]. The superfluid frac- 
tion /s vanishes much faster than the visibility V and 
the number fluctuations (Xj. For values of V/J where the 
superfluid fraction has practically vanished the visibil- 
ity is still larger than 0.7. Thus neither fringe visibility 
nor number fluctuations are a suitable indicator for the 
superfluid properties and the Mott transition in lattice 
systems. 

Quasi- Momentum Distribution. The interference 
pattern after ballistic expansion is closely related to 
the quasi-momentum distribution of the Bose gas in 
the lattice. Formally, the connection is revealed by 
constructing an expression for the occupation numbers 
for the Bloch states of the lowest band. We can use the 
relation between localized Wannier functions w{x — ^i) 
and Bloch functions ipq{x) to define a creation operator 
for a boson in Bloch state with quasi-momentum q 



4 = 77 E ^1' (13) 

i=l 

where £,i is the coordinate of the center of the zth lat- 
tice site. This relation is identical to the definition of the 
amplitude operator A{y) in (|l0| ) if we identify the phase 
(t^iiy) with q^i or the phase difference Scj) with qa, where 
a — Ci+i —£,1 is the lattice spacing. The occupation num- 
bers hq for the Bloch states with quasi-momenta q are, 
therefore, directly related to the intensity ( [TTI ) through 

= (*o|ctcJ*o)=W = 'Za). (14) 

Notice that in a finite system of length L the quasi- 
momentum q has discrete values which are integer mul- 
tiples of 2tt/L. The values of S(p = qa for these allowed 
quasi-momenta are marked by gray arrows in Fig. |^. 

Because of this intimate relation the interference 
pattern provides complete information on the quasi- 



momentum distribution of the trapped system. The in- 
tensity of the principal interference peak is proportional 
to the occupation number of the q = Bloch state, i.e., it 
describes the number of particles in the condensate. The 
washing out of the interference peaks with increasing in- 
teraction strength is linked to the successive redistribu- 
tion of the population from the condensate state with 
g = to states of higher quasi-momenta. In the limit of 
large V/J the intensity distribution is flat, i.e., all quasi- 
momentum states of the lowest band are occupied uni- 
formly. On the basis of this one-to-one correspondence 
between interference pattern and quasi-momentum dis- 
tribution we can reinterpret the visibility V of the in- 
terference fringes as measure for the uniformity of the 
quasi-momentum distribution. Vanishing visibility cor- 
responds to a completely uniform quasi-momentum dis- 
tribution, whereas visibility V = 1 means that at least 
one quasi-momentum state is unoccupied. 

Conclusions. We have shown that the matter-wave 
interference pattern observed experimentally contains all 
the information on the quasi-momentum distribution of 
the lattice system but no direct information on the su- 
perfiuid fraction. The behavior of the superfluid fraction 
shown in Fig. |l| depends strongly on the properties of 
the excitation spectrum, which enters through the sec- 
ond order contribution /s^^-*. The importance of this sec- 
ond order term shows that one cannot probe superfluidity 
through quantities which are only sensitive to the ground 
state of the system (like number fluctuation, condensate 
fraction, coherence properties, etc.). One has to devise 
an experimental scheme that measures superfluidity di- 
rectly. The formal deflnition of superfluidity gives a hint 
how to accomplish this. As mentioned earlier there are 
several methods to create the phase factor appearing in 
He experimentally, e.g. by accelerating the lattice. By 
observing the resulting flow behavior after release from 
the lattice one should be able to distinguish superfluid 
and non-superfluid components and determine the su- 
perfluid fraction. 
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